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RESEARCH MEMORANDUM 


INVESTIGATION AT MACH NUMBERS OF 1.62, 1.93, AND 2.41 OF 
THE EFFECT OF OSCILLATION AMPLITUDE ON THE DAMPING 
IN PITCH OF DELTA-WING—BODY COMBINATIONS 


By Arthur Henderson, Jr. 
SUMMARY 


In order to assess the validity of the linear theory, which accounts 
for neither viscous nor amplitude effects, tests were made of four delta- 
wing and body combinations at Mach numbers of 1.62, 1.95, and 2.41, and 


covering a Reynolds number range of O.H} x 106 to 1.69 x 106 to determine 
the effect of oscillation amplitude on the damping in pitch. 


The body in each case was the same; the wings had 259, 30°, 35°, 
and 45° semiapex angles and were flat plates with beveled edges. 


The results show that in some cases there is a small variation of 
damping with amplitude which is usually one of increasing damping with 
decreasing amplitude. А comparison of theory with experiment indicates 
that, in general, the damping in pitch of delta-wing—body combinations 
may be satisfactorily approximated by theory. 


INTRODUCTION 


The theoretical prediction of the damping in pitch of wings is 
limited by the restrictions inherent in the linear theory which, by its 
nature, is unable to show any effect of viscosity or variation of the 
damping-in-pitch parameter Сш + Сш, with amplitude. The responsibility 


of determining where, if at all, the existing theory is applicable, lies 
with the experimenter. There are two basic experimental techniques 
avallable by which the damping-in-pitch parameter can be determined, by 
the free-oscillation technique апа the forced-oscillation technique. 
References 1 and 2 present data obtained from free oscillations at both 
subsonic and supersonic speeds and reference 5 presents data obtained 
from forced oscillations at subsonic speeds. Reference 3 shows the 
effect of oscillation amplitude and frequency on damping in pitch. 
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The present paper presents the effect of oscillation amplitude on 
damping in pitch &t three supersonic Mach numbers, 1.62, 1.95, and 2.h1, 
for three wing-body combinations. А fourth configuration, the data for 
which were not reducible, is discussed. 


SYMBOIS 
с root chord (including part of wing submerged in the body) 
[d mean aerodynamic chord, 2/3c 
Co distance from wing apex to point of rotation 
D RETE moment per time rate of change of angle of attack 
1 moment of inertia 
K restoring moment per unit angle of attack 
1 distance from model center of gravity to point of rotation 
M Mach number or pitching moment 
m mass of model 
q angular velocity of pitch 
R Reynolds number 
S wing area (including part of wing submerged in the body) 
t time 
V velocity 
a angle of attack 
Ag amplitude 
В = М2 - 1 
Е wing semiapex angle 
р density 


д) frequency, radians / вес 
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Е! 


reduced frequency, ac/2V 
Cn pitching-moment coefficient, моиве 
с Sm 
Tag 
Cy, = n 
352 
2V 
ӘС 
. = oc 
"n = “Gg 
2ү 
Subscripts : 
1 wind off 
а wind on minus wind off 
cg center of gravity 


А dot above а symbol denotes differentiation with respect to time. 


APPARATUS AND TESTS 


Wind Tunnel 


All tests were made in the Langley 9-inch supersonic tunnel. The 
tunnel is of the continuous-operation complete-return type in which the 
stream pressure, temperature, and humidity can be controlled. Within 
the stagnation chamber ahead.of the first minimum are located 11 fine- 
mesh turbulence-damping screens. Тһе Mach number is varied by inter- 
changing nozzle blocks which form test sections approximately 9 inches 
square. The air was dried sufficiently for each test so that the con- 
densation effects in the test section were negligible. 


Models 


The models tested in the present investigation consisted of four 
delta-wing and slender-body combinations. The body in each case was the 
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same (fineness ratio, 11.2) апа consisted of а parabolic nose section, 

& circular cylindrical section along the wing-body juncture, and а boat- 
tailed afterbody. The wings had semiapex angles of 259, 30°, 35°, 

and 450, and the ratios of body diameter to wing span were 0.147, 0.132, 
0.120, and O.100, respectively. 


Figure 1 presents the models tested. Figure 1(а) gives the body 
dimensions and shows the relative wing positions and figure 1(Ъ) is a 
photograph of the component parts. The models were so constructed that 
each configuration consisted of four parts: (1) а nose апа hollow barrel 
with internally threaded end; (2) an internally and externally threaded 
collar which screwed into the hollow-barrel end; (5) а sleeve which slid 
over the hollow barrel and to which were soldered the wings to be tested 
and two smaller flat-plate yaw-stabilizing fins; and (4) a hollow shaft 
the threaded end of which screwed into the internal threads of the collar 
and the flat-head end of which was beveled to fair in with the rear of 
the boattailed sleeve. 


Model Support System 


In figure 2 is shown & model mounted in & mock-up of the test sec- 
tion. The support system consisted of & crossed flexural pivot assembly 
within the model and а 1/16-inch-diameter piano wire. One base of the 
flexural pivot (see fig. 5) made а push fit with the hollow-body barrel 
to which it was pinned when in place. The other base of the pivot was 
smaller in diameter than the inside diameter of the barrel, allowing 
this end to oscillate with respect to the body. The piano wire was 
fastened securely to one tunnel wall, passed through the free end of the 
pivot by means of oblong slots in the sides of the body, апа was secured 
by а yoke with а screw adjustment to а strain-gage beam on the outside 
of the other wall of the tunnel. (See fig. 2.) The screw adjustment 
and strain-gage beam allowed the tension on the wire to be held constant 
at approximately 250 pounds throughout all tests. Two set screws in the 
flat end of the free flexural pivot base clamped the base to the wire 
so the model was free to oscillate with respect to the wire. 


Mounted on the side wall of the tunnel behind the model was a spring- 
loaded mechanism through which the positioning of a guided wire could be 
controlled. The model was given its initial displacement by inserting 
the end of the guided wire into & small hole in the base of the model 
and retracting the wire any desired distance. Triggering the mechanism 
Successively released the wire from the body and removed the guide from 
the wake of the body. These two actions were practically simultaneous. 


Recording of Data 


A line source of light from а GE-855 bulb was focused on а 1/16-inch- 
diameter silvered mirror which was imbedded in the surface of the body 
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directly over the point of rotation. The line from the light source was 
then reflected onto а moving film strip through а horizontal slit in a 
20-foot high-speed film drum. Since the film speed was not constant, 
the pulses from a GE-AR8 argon timing light were recorded on the film 
strip along with the model oscillation record. A typical oscillation 
record is shown in figure 4. 


Displacements on the film strip were calibrated with measured angular 
displacements of the model. 
EXPERIMENT 


Scope of Tests 


The present investigation was conducted at Mach numbers of 1.62, 
1.95, and 2.41 and covered a Reynolds number range of O.l x 106 to 


1.69 x 106 based on the mean aerodynamic chord of the wing (including 
the portion submerged in the body). 


Every model was tested at two or more Reynolds numbers at each Mach 
number. In order to check repeatability, two or more complete sets of 
data were taken and all are presented. All models were oscillated about 
ОО angle of attack. 

Reduction of Data 

The single degree of freedom, free-oscillation technique was used 

to determine the damping-in-pitch characteristics of the models of the 


present investigation. With viscous damping and а nonvariable spring 
constant, this motion is described by the following equation: 


Ta + Da + Ka = 0 (1) 
The solution to equation (1) is 
- Dx 
2. el D 2 
с (со) 408 (сов ай + PTa sin ot) ( ) 


where 


2 
бу = К (D 
(а) 
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а = age PI 2 


In the ideal case the term = (which is a measure of the damping) is 


determined by plotting a of the envelope against t ona semilog scale. 
This plot is a straight line, the slope of whichis - = 


In the present investigation, since the oscillation-record envelopes 
could not be represented by equation (3) throughout the whole oscilla- 
tion, it was assumed that the oscillation record could be represented by 
& series of segments, each segment of which could be represented by equa- 

D 


tion (3) with the appropriate value of aT" These segments may be made 
as small as desired. Differentiating equation (5) with respect to time 
and letting time approach zero for each segment gives 


En өнд e) 


Thus - ів the reciprocal of the number of time units lying between 


the time ordinate of & point on the envelope and the intercept on the 
a= О axis of the tangent to that point. The term = is then deter- 


mined as а function of amplitude. 


To make use of the previously mentioned method of reducing the data, 
the following procedure was employed. The peaks of the oscillation 
record were read on а telereader and plotted. The amplitudes were plot- 
ted to 5 times actual size and the time scale was adjusted to facilitate 
plotting. Curves were then faired through the upper and lower peaks 
and the а = О line was determined. At arbitrarily selected points 
along the upper curve, tangents were drawn to the curve and extended 
until they intersected the а = О line. The same thing was done at 
corresponding points on the lower curve. Тһе reciprocal of the number 
of time units between the time ordinate on the curve and the intercept 


on the а = О line was taken to be T for the point. 
In the present paper, both the damping-in-pitch parameter Cmq + Спа, 
and the pitching-moment derivative Сы, are found. Inasmuch as, with 


the wind on, equation (1) is 
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‚т | 
14 + (Da #D4)& + (Ka + Kija = 0 (5) 
it is necessary to determine I, Di, and K,, in order to be able to 


separate the aerodynamic terms from the total result, and to evaluate D 


D 
fr — o 
om DI 


The moment of inertia at the center of gravity Icg was measured 
according to the bifilar pendulum equation 


= шеа = г) 
Ieg x (6) 
where 
W weight of configuration 
a) frequency in radians per second 
а, в, г symbols defined in figure 5 


The value of I for the test was then determined by the equation 


Т = Ieg + mî 


Equation (6) is derived on the assumption of small angles; that is, sines 
and tangents of angles are replaced by the angles themselves. Therefore, 
s Should be long compared with d, and the value of r should be as 
close to 0.5 as practicable in order to prevent a coupling of gravity 
and torsional-pendulum effects. 


K; was determined by the equation 
Кі = Там? (7) 
where аң is the wind-off frequency. 


D; was determined from the equation 


Dy = 212) (8) 
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where 2 is the wind-off damping factor. 
1 


It might be mentioned here that the seguence of tests for each run 
was as follows: one wind-off test, two tests at one Reynolds number, 
two tests at the second Reynolds number, and a second wind-off test. 


Since Кі and TD, were found to be essentially invariant with 


pressure from 1 to 0.1 atmospheres and there was no significant change 
in Пі with amplitude, the Кі and Dj values used were the mean of 


811 wind-off tests. In equation (5), 


Од = ӘМ + eM 
да да 
- lyss Cm + две Sm (9) 
2 да, 2 да 
Inasmuch as in stability notation Спа апа . аге defined аз Om 
ac 
d 
and Cn , respectively, 
g 9E 
ev 
Са + Ong, = № (10) 
pSVc 
Similarl = Cm 
у, since Cy, = 5 
Ka 
ы 11 
Cir: pe (11) 


NACA RM 153H25 (Som. 7) р B 


The value of Ка was determined by subtracting К: from the total wind- 
on 8pring constant obtained from 


Precision of Data 


There are two main sources of error involved in thís method, one 
source is directly associated with the flexural pivot and the other is 
inherent in the data reduction method. 


With regard to the flexural pivot, it was originally intended to 
use several different flexure plates of varying thicknesses in order to 
attain a fairly wide range of frequencies for the damping tests. As the 
thickness of the plates increased, however, not only did Кі and р; 


increase, they also became increasingly erratic. This condition was a 
direct result of size limitations. The dimensions of the flexural-pivot 
assembly (fig. 3) were such that 0-80 screws were the largest that could 
be used to hold the flexure plate rigidly against the flexure-plate base. 
This size was not sufficient to prevent some friction from taking place. 
To lessen this effect, only the thinnest flexure plates (0.008 inch thick) 
were used. The characteristics of these plates, as determined from all 
wind-off tests were: 


(2) = 0.042 + 0.007 


21/1 
which was 11.2 percent of the smallest = value found from all wind-on 
tests, and 

K4 = 0.18 + 0.002 
which was 11.3 percent of the smallest K value for the test. 


With regard to the errors involved in the data reduction, it is nec- 
essary to determine the distance from a point on & curve to the inter- 
section of the tangent to the curve at that point with a straight line. 

As the angle between the tangent and the straight line decreases, the 
error involved in the graphical construction increases. Thus, the prob- 
able error involved varies with angle of attack, the error being least 

at the higher angles of attack and greatest at the lower angles of attack. 
In view of these conditions, it is difficult to determine some particular 
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number which would represent the overall precision of the data. Since 
the previously mentioned source of error probably overshadows all other 
probable errors it is believed that it should merely be pointed out that 
Спа + Cm; is more accurate at the higher angles of attack than at the 


lower ones, and that апу peculiar behavior of Ста + . at the low 
angles of attack should be viewed with suspicion. 


For the sake of completeness, the following additional uncertainties 
are pointed out: 


(а) Although, as indicated in reference 5, Dj may vary with fre- 
quency, no information on the possible variation of Р; with frequency 
was obtained for this series of tests. 


(b) The value of D, is determined under no-load conditions and 


this quantity is then subtracted from the total wind-on damping when the 
Spring mechanism is subjected to oscillating lift and drag forces. It 
is believed that the number used to represent Dj is at least of the 


correct order of magnitude; however, actual load conditions were not 
simulated. 


(c) The support wire, although under 250 pounds of tension, is none- 
theless flexible and therefore subject to deflection under load.  Esti- 
mates of the drag of the 209 delta-wing—body combination at R = 1.69 x 106 
give а zero lift drag of about 5 pounds plus an oscillating drag force 
of from О to 1 pound for a range of amplitude of 0° to 59. A dead weight 
of 5 pounds suspended from the center of the support wire gave a deflec- 
tion in the wire of 0.052 inch. Assuming а linear displacement with load, 
as the model oscillates, the center of the wire will move rearward period- 
ically a maximum of about 0.01 inch. It may therefore be said that the 
center of rotation is known to within less than 0.01 inch, this amount 
decreasing as the amplitude decreases. 


RESULTS 


Test results on the 45° wing were irreducible and hence are not pre- 
sented. There are two related reasons for this irreducibleness. 


The first reason is associated with the airfoil shapes of this series 
of tests, which were flat plates with beveled edges. А simple two- 
dimensional analysis shows that this type of section, if there is any 
tendency toward flow separation, has two positions of zero lift, ta, 
where a is a function of the pressure-distribution change associated 
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with the separation. Figure 6 shows part of the oscillation record for 


the 45° wing at М = 1.95 and R= O.5 x 10° from which it can be 
seen that, although the wing initially oscillated about 0° angle of 
attack, after a few cycles it suddenly began to oscillate about some 
angle of attack other than zero. 


с 
The second reason is that for the 150 wing = 0.94; that is, the 


point of rotation is located at 0.О4с ahead of the theoretical center of 
pressure for a flat plate. With such a small moment arm, the effect of 
any shift of the center of pressure will be greatly magnified and will 
show up ав a variation in frequency, which is a measure of K4 + Ка. 


Figure 7 is a plot of frequency in cycles per second against time, as 
taken from figure 6. Superimposed on this plot is the natural frequency 
of the model with the wind off. The difference between these two curves 
ів а measure of Cm," ТЕ can be seen that while the model was oscillating 


about 0° angle of attack, Cm, decreased with time from positive to neg- 


ative values; that is, as a decreased, the center of pressure moved 
from behind to ahead of the point of rotation. Then, when the mean oscil- 
lation angle changed from zero to its new value, Сто, became discontin- 


‘uously positive and increased with time, or there was a discontinuous 
rearward shift in the center of pressure, after which it continued to 
move rearward. Thus, the separation coupled with the small moment arm 
of the 45° wing combined to make the data irreducible. For the tests of 
the 45° wing at higher Reynolds numbers the effects mentioned here were 
less pronounced, but the tendencies were the same. 


Although the same wing section was used for the other wings of the 
investigation, it is believed that the separation tendency noted on the 
45° wing did not interfere with data reduction for the other wings because 
for the 250, 30°, and 35° wings the mean oscillation angle did not change 
from 0°. 


Figures 8, 9, and 10 present the variation of Cmq + Cms with ampli- 
tude for the wings of 25°, 30°, and 55° semiapex angles, respectively. 
Figures 11 to 15 present the corresponding plots of Cm, against amplitude. 


The solid lines on figures 8 to 10 represent the theory of reference 4, 
whereas, the solid lines on figures 11 to 15 represent the theoretical 
prediction obtained by combining the method of reference н with that of 
reference 5. 


Each plot of Cm, + Cmg, with the exception of figure 10(b), has at 


least four sets of points representing two or more different tests. For 
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each test the unflagged symbols represent values obtained from the top 
oscillation curve and the flagged symbols répresent values from the 
bottom curve. Where an unflagged symbol occurs without а corresponding 
flegged one, the two points coincided. 


Figures 8(a) and 9(a) contain more sets of points than the other 
figures because some unusual results of the original data were questioned. 
Additional tests at а higher initial amplitude than originally employed 
were made which showed that in а certain frequency range and above an 
amplitude of about 1.59 the support system interfered with the model. 

This interference was evidenced by the appearance of beats in the oscil- 
lation record in the higher amplitude range. The beats were a clear 
indication that some external source was feeding periodic pulses of energy 
back into the system. This source could only be the support system. 
Figures 8(a) and 9(a), at the lowest Reynolds number, fell within this 
frequency range. Therefore data above about dp = 1.59 are not presented 


for these figures. 


On the left-hand side of figure 8(a) are presented four sets of 
data; three sets at R = 1.31 x 106 and one at R = 1.67 x 10°. Increasing 
R by 27 percent produced а noticeable change in the damping in the higher 
amplitude range апа caused а greater variation of damping with amplitude. 
The three sets of data on the right-hand side of figure 8(a) check fairly 
well except in the very low amplitude range, where, as previously pointed 
out, the accuracy is doubtful. 


The data of figures 8(c) and 11(с) are questionable. During the test 
at В = 0.96 x 10” the model was fluttering slightly in yaw. As the 
Reynolds number was decreased both the yawing amplitude and frequency 
increased. At В = 1.0h x 100 the flutter ceased abruptly. The next 


two tests were made at a nominal Reynolds number of 1.09 x 106; the 
first of these tests (diamond-shaped symbol) was actually at В = 1.07 x 106 


and the second (triangular symbol), at 1.10 x 106. It is difficult to 
conceive of such а small change in Reynolds number causing such а large 
difference in the two curves shown. Although these two tests were made 
close to the yaw stability boundary, the fact that the lower Reynolds 
number test, which was in the flutter-in-yaw regime, did not exhibit 
any drastic changes in damping suggests that the changes noted are not 
due to any possible coupling of lateral апа longitudinal effects. 


This particular test took place in about the middle of the series 
and inasmuch as all succeeding tests were satisfactory, equipment trouble 
is eliminated. Since the wing leading edge in this case was transonic 
(B tan e = 1.02h), figure 8(c) may be an indication that when the Mach 
lines аге close to the leading edge, а random variation of the damping 
can be expected. 
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Figure 9la) on the left consists of three sets of data: two at 
R= 1.17 x 10”, one with and one without transition strip, and one at 


В = 1.69 x 106 without transition strip. In the upper amplitude region 
&t В = 1.69 x 106, the damping is slightly lower than for the other two 
tests, increasing with decreasing amplitude until it joins, and tbere- 
after follows, the test at В = 1.17 x 100 with transition strip. The 
damping for R = 1.17 x 106 without transition strip develops а hump 

in the vicinity of оо = 0.750. Inasmuch as this hump approaches the 


level of the damping for R = 0.75 x 106 (on the right-hand side of 

figure 8(a)) without transition strip, it would appear that the hump is 
caused by separation effects in the lower amplitude range. It is note- 
worthy that the test at R = 0.75 x 106 with transition strip gives а damping 
which is of the same order of magnitude as that for the higher Reynolds 
number tests. 


И 


| 


А general examination of figures 8 to 10 indicates that within the 
range of these tests, Cmq * Сте, may be satisfactorily approximated by 


theory. In some cases, the experimental results indicate a slight vari- 
ation of Сша + Сте, with amplitude. This variation, when it occurs, is 


generally one of increasing damping with decreasing amplitude. 


From figures 1l to 15 it can be seen that the accuracy of the theo- 
retical prediction of Cm, increases as p tan е increases until, when 
the leading edge is well supersonic, the agreement between theory and 
experiment is excellent. 


CONCLUDING REMARKS 


An investigation has been conducted at Mach numbers of 1.62, 1.93, 
and 2.41 to determine the damping in pitch of four delta-wing and body 
combinations Рог a Reynolds number range of 0.11 x 106 to 1.69 x 106 
based on the mean aerodynamic chord of the wings. The semiapex angles 
of the wings were 25°, 30°, 35°, and 45°. Тһе results are presented as 
a function of amplitude. Within the scope of these tests the following 
conclusions are indicated: 
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Тһе damping-in-pitch parameter Сша + Cm шау be satisfactorily 
approximated by the theory. 


Where & variation of Ста + Cg; with amplitude occurs, it is usually 
small and the damping increases with decreasing amplitude. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., August 18, 1953. 
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Figure 1.- Test models. 
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(b) Component parts of models. 1-80292 
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Figure 1.- Concluded. 
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Figure 2.- Model mounted in & mock-up of the test section. 
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Figure 5.- Two views of flexural-pivot assembly. 
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Figure 4.- Typical oscillation record; е = 25°; Ме 1.62; В = 1.51 X 10°; 
initial цо = 2.70°. 
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Figure 5.- Bifilar pendulum setup. 
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Figure 6.- Section of oscillation record for є = 45°; M = Es 95; 
В = 0.54} x 10 
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Figure T.- Variation of frequency with time for 45° wing. М = 1.95; 
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(a) М = 1.62; В tan є = 0.59. 


(b) М = 1.95; В tan є = 0.770. 
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Figure 8.- Variation of Cm, + Cm, with amplitude for в = 259; 
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(а) M=1.62; 8 tan е 


(b) М= 1.95; В tan е = 0.953. 


(c) M= 2.41; В tan е = 1.266. 
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Figure 9.- Variation of Cm, + Cm; with amplitude for є = 30; = = 0.80. 


(Flagged symbols refer to bottom oscillation curve.) 
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Figure 10.- Variation of Cm, + Cm; with amplitude for Е = 559; = 0.85. 


(Flagged symbols refer to pottom oscillation curve.) 
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Figure 11.- Variation of Cm, with amplitude for 
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Figure 15.- Variation of Сто with amplitude for 
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